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Abstract 

The Dirac operator describing the couphng of continuum quark fields to SU{2) 
center vortex world-surfaces composed of elementary squares on a hypercubic lattice 
is constructed. It is used to evaluate the quenched Dirac spectral density in the ran- 
dom vortex world-surface model, which previously has been shown to quantitatively 
reproduce both the confinement properties and the topological susceptibility of SU (2) 
Yang-Mills theory. Under certain conditions on the modeling of the vortex gauge field, 
a behavior of the quenched chiral condensate as a function of temperature is obtained 
which is consistent with measurements in SU (2) lattice Yang-Mills theory. 
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1 Introduction 

The random vortex surface model is an infrared effective theory which aims to describe 

the relevant physical fluctuations determining the low-energy nonperturbative sector of the 
strong interaction. It is based on collective magnetic vortex degrees of freedom ||^-[^ which 
are represented by closed two-dimensional world- surfaces in four-dimensional (Euclidean) 
space-time. The chromomagnetic flux carried by the vortices is quantized according to the 
center of the gauge group. The vortex world-surfaces are treated as random surfaces, an 
ensemble of which in practice is generated using Monte Carlo methods on a hypercubic 
lattice; the surfaces are composed of elementary squares (plaquettes) on that lattice. The 
spacing of the lattice is a fixed physical quantity, related to an intrinsic thickness of the vortex 
fluxes, and represents the ultraviolet cutoff inherent in any infrared effective framework. This 
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model can be adjusted such as to quantitatively reproduce the confinement properties of 
SU{2) lattice gauge theory including the finite-temperature transition to the deconfined 
phase. It simultaneously predicts a quantitatively correct topological susceptibility, again as 
a function of temperature 0. The properties of the model vortices closely parallel the ones 
obtained for vortex configurations extracted from the full SU{2) lattice Yang-Mills ensemble 
using an appropriate gauge fixing and projection procedure 0,0, cf. ||TT[|-|[T^. 



The aim of the present work is to investigate whether also the third central nonperturba- 
tive phenomenon determining low-energy strong interaction physics, namely the spontaneous 
breaking of chiral symmetry, can be correctly described within the random vortex surface 
model. For this purpose, a number of technical developments are carried out. The coupling 
of the gluonic vortex degrees of freedom to (continuum) quark fields is governed by the Dirac 
operator; in order to construct this operator, an explicit representation of vortex configura- 
tions in terms of continuum gauge fields is needed. Having specified such a representation, 
matrix elements of the Dirac operator can be calculated in a truncated infrared quark ba- 
sis. Via the (ensemble-averaged) spectral distribution of the resulting Dirac matrix, which 
is obtained numerically, the (quenched) chiral condensate is extracted; this condensate can 
be used as an order parameter for spontaneous chiral symmetry breaking. Under certain 
conditions on the explicit vortex gauge field modeling, the chiral condensate as a function of 
temperature behaves in a manner consistent with SU (2) lattice gauge theory. 

To date, several explorations of spontaneous chiral symmetry breaking within the vortex 
picture have been presented. One of these lines of investigation |T5[-|T^ deals with the 



vortex configurations extracted from the full SU{2) lattice Yang-Mills ensemble by gauge 
fixing and projection Removing these vortices from SU{2) lattice configurations leads 



to the restoration of chiral symmetry |15|; on the other hand, the extracted vortices by 



themselves break chiral symmetry at zero temperature |T6[. In the chiral limit, a divergence 



of the chiral condensate is observed which is reminiscent of the behavior detected in full 



SU{2) Yang- Mills theory using domain- wall fermions [|T8|; such a divergence will also be 
observed in the Dirac spectra obtained in the present work, cf. section |^. Crossing into the 
deconfined regime, chiral symmetry remains broken, to an extent which differs qualitatively 
| |T7| from the breaking found in the aforementioned domain-wall fermion calculations. In view 
of this result, the significance of the chiral condensate induced by vortices and its relation 
to the condensate obtained in full Yang-Mills theory remained in doubt. In the present 
work, the behavior of the chiral condensate in the high-temperature regime will turn out 
to depend strongly on the modeling of the vortex gauge field. Only a sufficiently smooth 
definition of the gauge field yields a behavior of the chiral condensate consistent with full 
Yang-Mills theory. Possibly the vortex extraction procedure used in endows the vortices 
with excessive spurious ultraviolet disorder which leads to the unnatural result for the chiral 
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condensate in the deconfined phase. 

On the other hand, quite recently an effective vortex model (in the guise of a Z{2) 
gauge theory in the strong coupling limit) was considered which to a certain extent 



follows the logic of the random vortex surface model investigated here. At the expense 
of a more limited scope of applicability (the ultraviolet cutoff is e.g. too low to permit a 
description of the deconfinement phase transition), the dynamics is simplified to the point 
where analytical estimates of the chiral condensate become possible. These estimates appear 
to be close to phenomenologically accepted values for the unquenched, dynamical quark case. 
However, the extreme limit considered seems rather artificial in that the quenched and the 
dynamical quark cases yield comparable chiral condensates. Presumably, since one possibility 



of arriving at estimates of the type used in ||T9[ appears to lie in invoking the large-A^^c 
expansion {Nc denoting the number of colors), this should not come as a surprise; the large- iVc 
expansion is by construction quenched to leading order, dynamical quark effects only entering 
subdominant terms. By contrast, the quenched chiral condensate obtained in lattice Yang- 
Mills theory |^ exceeds the one expected phenomenologically in the presence of dynamical 
quarks by about an order of magnitude; clearly, the quark determinant weighting of the 
dynamical quark ensemble must be viewed as one of the dominant physical effects in order 
to reconcile the quenched approximation with the dynamical quark case. In view of this, 
the meaning of the chiral condensate discussed in remains unclear, as does the question 
whether the quoted values in fact support or weaken the basis for that model. 



2 General remarks on the formulation of a vortex gauge 
field model 

In the applications of the random vortex surface model hitherto considered it was not 

necessary to explicitly construct a gauge field associated the vortex configurations arising in 
the model. To evaluate confinement properties, encoded in Wilson loops fl]], it is sufficient 
to specify the space-time location of the vortex world- surfaces; to evaluate the topological 
charge 0, it is sufficient to specify the vortex field strength, which is localized on those 
surfaces. By contrast, construction of the Dirac operator implies specifying the gauge field 
itself. The freedom inherent in this specification entails ambiguities^ due to the fact that the 
concept of an infrared effective model of gauge and quark fields to a certain extent clashes 
with gauge invariance; after all, the distinction whether a given configuration is infrared 

^Further to the ambiguity discussed in the main text, an additional one in principle can arise from the 
fact that, given a particular (vortex) field strength, in the full space of nonabelian gauge fields, this field 
strength may be encoded by several gauge-nonequivalent gauge fields |^ . The latter however will in general 
yield different Dirac spectra. 
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or not is a gauge-dependent one. In an infrared effective model of gaiige and quark fields, 
manifest gauge invariance can only be sustained with respect to gauge transformations U 
wfiich are themselves infrared. Specifically, when using a truncated basis of quark wave 
functions (as will be done in this investigation), invariance under a change of the gauge in 
which the gauge field ensemble is given is only guaranteed to the extent that, for a quark 
wave function q represented in the aforementioned basis, also Uq can be represented to good 
accuracy in that basis. Conceptually, this is not an obstacle to model-building; it just implies 
that a full specification of the model must include a fixed choice of gauge, and that predictive 
power is curtailed to the extent that observables depend on the precise gauge field modeling. 
In general, some observables have to be used to fix the model before additional ones can be 
predicted. 

Ideally, the most consistent gauge field description of a given vortex field strength ensemble 
would be one in which the gauge field is itself as smooth, i.e. infrared, as possible, such 
as the Landau gauge. In practice, the choice is constrained by considerations of technical 
manageability, and also by certain idealizations of the vortex configurations which already 
occur on the level of the field strengths, cf. below. In order to obtain information about the 
dependence of observables on the construction of the gauge field, in this work different choices 
will be explored, which, configuration by configuration, describe the same vortex world- 
surfaces, but are related by gauge transformations which by no means are infrared. It will 
turn out that chiral symmetry breaking in the confined phase is quite robust as the modeling 
of the gauge field is changed; on the other hand, the Dirac spectrum in the deconfined, 
high temperature phase exhibits a substantial dependence on the gauge field modeling - 
gratifyingly, the most consistent (smoothest) construction turns out to best reproduce the 
results of full Yang-Mills measurements. 

Before proceeding to present the construction of the gauge field in detail, a few remarks 
are in order concerning the idealizations of the vortex field strengths mentioned further 
above, which were already introduced in previous considerations of the random vortex surface 
model. The vortices hitherto used to evaluate observables within the model were idealized 
in two ways: For one, vortex fiuxes were taken to be infinitely thin, i.e. the associated 
world-surfaces are truly two-dimensional, with no transverse thickness in the remaining two 
space-time directions. Note that this idealization is gauge-invariant; it is a statement about 
the structure of the modulus of the field strength trF^. On the other hand, as a vector in 
color space, the field strength F on vortex surfaces was taken to point exclusively into the 
positive or the negative 3-direction; generic vortex surfaces are non-oriented, i.e. consist of 
patches of alternating orientation. Note that the color direction of F can be changed by 
gauge transformations; thus, the configurations hitherto used can be viewed as having been 
specified in a gauge of the maximally Abelian type. 
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It should be emphasized that these are the properties of the field strengths used to evaluate 
observables within the model. By contrast, in the generation of the vortex surface ensemble, 
a physical thickness of the vortices is indeed implicit : A finite action density is associated 
with vortex surface area and curvature; infinitely thin vortex fiuxes of course formally would 
carry divergent trF^. Also, vortex surfaces cannot be packed arbitrarily densely within the 
model; instead, there is a fixed ultraviolet cutoff, interpreted as a consequence of the vortex 
thickness. The construction and interpretation of the vortex dynamics is discussed at length 
ini. 

The idealized thin vortex configurations described above are, however, entirely sufficient 
for an unambiguous evaluation of the observables considered to date. Neither the asymptotic 



string tension |l|],p2[ nor the topological charge depend on any assumption about 



the transverse structure of the vortex fiuxes, i.e. whether they are thickened or not. Also, 
both observables are virtually independent of assumptions about the color orientation of the 
vortex field strength. For the topological charge, this is somewhat less obvious than for the 
Wilson loop; this issue is discussed in detail in 0. 

On the other hand, in the case of the Dirac operator spectrum to be discussed here, 
the issue is a priori not as clear-cut; it is therefore necessary to be specific about how the 
aforementioned idealizations will be dealt with. In the following treatment, the idealization of 
thinness of vortex fiux will be maintained; a sufficiently simple and manageable construction 
of manifestly thick vortex fiuxes did not become apparent in the course of this investigation^ 
Instead, a model assumption will be used which is familiar e.g. from the Nambu-Jona-Lasinio 
(NJL) quark model. Namely, as already hinted at above, the quark modes will be subject to 
a fixed ultraviolet cutoff, the same one which is implicit in the model vortex dynamics via the 
spacing of the lattice on which the vortex surfaces are defined . In the NJL model, the quark 
current-current interaction formally is point-like; however, the ultraviolet cutoff on the quark 
modes effectively smears out the four-quark vertex into a nonlocal interaction. In complete 
analogy, within the present treatment, the interaction between a quark and a vortex formally 
takes place on an infinitely thin submanifold (to be specific, a three-dimensional volume, the 
two-dimensional boundary of which precisely gives the location of the vortex field strength, 
cf. section [3.11 ). However, the ultraviolet cutoff on the quark modes will effectively smear this 
out in the direction perpendicular to the volume (and, at its boundary, in the two directions 
perpendicular to the idealized thin vortex field strength). Thus, the thickening of the vortices 
is not carried out explicitly in the vortex gauge fields, but is effectively accomplished by 
smearing out the vortex-quark vertex via the ultraviolet cutoff on the quark fields. 

•^It should be noted that the thickness of realistic vortex fluxes (about 1 fm in diameter ||^,|9j) noticeably 
exceeds the typical distance separating neighboring vortices; a two-dimensional plane in space-time is, on the 
average, pierced by 1.8 vortices/fm^ at zero temperature Physical vortex fluxes overlap to a considerable 
extent. 
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Turning to the issue of the color orientation of the vortex field strength, a more sophisti- 
cated approach will be pursued, for the following reason. If one works with a globally defined 
field strength, pointing either into the positive or the negative 3-direction in color space, 
and one attempts to globally define an Abelian gauge field generating this field strength, 
then, in the (generic) case of nonorientable vortex surfaces, one necessarily introduces Dirac 
string world- surfaces into the description, emanating from the lines on the vortex surfaces at 
which the orientation switches. After all, a globally defined (Abelian) gauge field must satisfy 
continuity of (Abelian) magnetic flux. This would not constitute a problem if one worked 
with a complete quark basis, i.e. solved the Dirac equation exactly. Then, the quark wave 
functions would exhibit singularities along the Dirac strings which would serve to cancel any 
physical effect of these strings, i.e. the Dirac strings would be unobservable, as they should 
be. However, in any truncated calculation such as will be pursued here, the cancellation 
would not be perfect; instead, the Dirac strings effectively would act as additional physical 
magnetic fluxes (of double magnitude compared with vortex fluxes) with which the quarks 
can interact. Thus, more magnetic disorder would be effectively present than the model aims 
to describe. 

The way out of this dilemma lies in instead using a Wu-Yang construction of the gauge 



field ||2T|. In other words, the gauge field will be defined on local patches, which are then 



related to one another via transition functions on their overlaps; the transition functions 
in general will have to be non-Abelian. The patches should be chosen sufficiently small, 
such that, on any given patch, the vortex surfaces can be oriented; in this case, no Dirac 
strings arise in the associated Abelian gauge field, which is then well-defined on the whole 
patch. The reader should however be forewarned that constraints on the space-time form 
of the transition functions in generic vortex configurations will force the reintroduction of 
some spurious fluxes (via additional non-Abelian gauge field components) which are formally 
similar to Dirac strings, but with the crucial difference that they manifestly decouple from 
the smooth infrared quark fields. Thus, spurious fluxes are not completely eliminated by the 
construction presented below; rather, they are recast in an innocuous form. In practice, the 
Dirac equation will be solved using the finite element method, i.e. using a basis of localized 
quark wave functions. Thus, it is natural (and sufficient) to use the supports of the individual 
basis functions, i.e. the finite elements, as the space-time patches on which the gauge field is 
locally defined. Specifics on these finite elements follow further below. 
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3 Construction of the vortex gauge field 



3.1 Local structure of the gauge field 

A closed SU (2) vortex world-surface S along with its quantized chromomagnetic flux can, to 
a large extent, be characterized by the gauge- invariant property 

W[C] = (-l)^(^'^) (1) 

for any Wilson loop W[C] linked L{C, S) timesQ to the surface S. Note that the factor (—1) 
on the right hand side of eq. (|l|) corresponds to the only nontrivial center element of the 
SU{2) group, which will be the subject of investigation in the following. It is in this sense 
that vortex flux is quantized according to the center of the gauge group. For higher SU{N) 
groups, the (—1) on the right hand side of eq. (|ip is replaced by center elements of those 
groups. In general, there are several such center elements and, consequently, several types of 
quantized vortex flux. 

For the purpose of constructing a gauge field encoding the property (|^), it is useful to shift 
the emphasis of the description |^ , namely from the closed surfaces S to three-dimensional 



volumes S bounded by S, i.e. S = 9S. Then (P can be rewritten as 

W[C] = (-l)^(^'^) (2) 

in terms of the intersection number^ I{C, S) of C with S. 

The random vortex surface model studied here is formulated on a space-time lattice, on 
which vortex surfaces are composed of elementary squares; correspondingly, the volumes E 
are composed of elementary three-dimensional cubes. To specify a volume E, as yet without 
orientation, it is sufficient to associate each elementary three-dimensional cube K on the 
lattice with a value 4>{K) = ±1, (^{K) = 1 signifying that K is not part of E, and (piK) = —1 
signifying that it is. The question of how a volume E for a given vortex surface S can be 
obtained in practice will be discussed further below; for the time being, let it be assumed 
that such a volume has been found (obviously there is some freedom in choosing E, since it 

*For physical vortex fluxes, which are endowed with a certain thickness in the directions perpendicular to 
S, the contour C would have to circumscribe S at distances larger than this thickness in order to capture 
the entire vortex flux and not intersect it. As already discussed further above, in the present treatment, no 
explicit transverse vortex thickness will be introduced. Thus, such a qualification does not arise. 

^In the exponentiated form (^), it is irrelevant whether one counts intersections without reference to the 
relative orientation of E and C, or weighted with a sign according to the orientation. Below, when defining 
the gauge field itself, it will be necessary to be more specific, and the orientation will play a crucial role. 
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merely must satisfy (9S = S). The Wilson loop (Q) in this language become^ 

W[C]= n <PiK.) (3) 

where the product is taken over all cubesQ Ki intersected by C. 

This must now be translated into the gauge field A{K) associated with each cube K, i.e. 
a gauge field must be constructed which satisfies 

exp(tJ^dx^A^{K)^ =<j){K) (4) 

for a path P in /i-direction intersecting the three-dimensional cube K, assumed here to extend 
into all directions but the /i-direction. Gauge fields A^^^K) which reproduce this property 
can be given as follows. If K is the cube extending from the lattice site x into all (positive) 
directions except the /i-direction, then (with a denoting the lattice spacing and the third 
Pauli matrix encoding the color structure) 

A^{K) = n{K)7ra^6{x^ - x^) J] 0{xx - Xx)e{xx + a - xx) (5) 

i.e. A{K) has support only on K and points into the space-time direction perpendicular 
to K] furthermore, its strength depends on an integer n{K) which, at this point, is merely 
constrained to satisfy 

exp{iTm{K)) = (piK) , (6) 

i.e. n{K) must be even for (f){K) = 1 and odd for (f){K) = — 1. Apart from this, the integer 
n{K) still remains to be specified (which is where the orientation of S will attain relevance, 

^Each cube K corresponds to a link on the lattice dual to the one used in the present construction, 
that link perpendicularly piercing the cube K in question. This dual lattice is the lattice one would use to 
define standard lattice gauge theory; indeed, if one transfers the values (j){K) to the corresponding dual links, 
one has precisely constructed a Z{2) lattice gauge configuration with the vortex content given by S. The 
Wilson loop is then the product over all links it contains, as is manifest in eq. (|^). Nevertheless, it should 
be emphasized that the gauge fields defined in the following, designed such as to reproduce (||), contain 
considerably more information than the aforementioned Z{2) lattice configuration. This information, as will 



become clearer in section 3.2, is encoded in the manner in which different cubes Ki are sewn together at 
their shared boundaries such as to compose entire three-dimensional volumes E. These volumes contain 
unambiguous global topological information, reflected in the gauge fields constructed here; by contrast, in 
the standard lattice formulation, such topological information is, at best, very implicit. The gauge fields 
obtained in the present construction are continuum gauge fields, even if they do assume a decidedly cubistic 
form. 

^In the following, when more than one cube K is being referred to, this will be indicated by supplying an 
index label i, i.e. Ki, where it should be clear from the context which set the index i runs over. 
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along with the necessity of dividing space-time into patches, as already indicated in section 0). 
Finally, the color structure of the gauge field is given by the third Pauli matrix in (^). As 
already discussed in section 0, this is a specification which was already adopted previously 
0,p3| on the level of the vortex field strengths; it implies that vortex configurations at 



this stage are cast in a gauge of the maximally Abelian type. Another strong restriction of 
the gauge freedom is of course implied by the space-time form of (|^); the gauge field only 
has support on a three-dimensional submanifold made up of elementary three-dimensional 
cubes in the lattice^, as already indicated and discussed in section ||. In fact, apart from 
the integer n{K) to be specified in (^, the only freedom which remains is the choice of the 
three-dimensional volume S spanning the given vortex surface S. This freedom corresponds 
to the Z{2) part of the SU{2) gauge freedom of SU{2) Yang- Mills theory pf. A gain, the 
question of how E can be obtained for a given vortex surface 5* in practice is deferred for 
now; suffice it to state that more than one choice will be explored below. 



3.2 Global topology and space-time patches 

In a gauge field formulation of vortex configurations, one needs a more specific characteri- 
zation of vortices than is contained in the property ([l|); this formally manifests itself in the 
a priori freedom in the choice of the integer n{K) in (^. Within the random vortex surface 
model, the only physical degrees of freedom are elementary vortex fluxes whose gauge field 
description satisfies 

dx^Af, = ±na^ (7) 

for a line integral along an infinitesimal loop C circumscribing the vortex]^. 

This property will serve to fix the integer n{K) in (^) to a large extent. Eq. (|^) should be 
contrasted with the exponentiated form (|l|), which would allow any odd integer multiple on 
the right hand side[^ of (|^). It should be noted that the specification (0) was already used 
when calculating the topological charge of vortex configurations @,|2^; fluxes corresponding 

^This general type of gauge field was called "ideal vortex field" in |2^, with the slight difference that in 
p3| , the ideal vortex field had support only on E, i.e. on the subset of cubes Ki with 0(-fi'i) — —1. Below, the 
definitions will in fact be brought to match even more closely, by viewing those Ki which are associated with 
4>{Ki) — f , but which nevertheless support a nontrivial gauge field (such Ki are unavoidable if one wishes to 
exclude Dirac strings, cf. section ^ ), as two cubes Kii,Ki2 with ^{Kn) = 4>{Ki2) = — f which are parallel 
and slightly displaced from one another; i.e., these are viewed as additional parts of S which initially could 
not be resolved due to the coarse-grained character of the lattice description. 

^Of course, for physical thick vortices, the relevant loops would be ones of a size corresponding to the 
vortex thickness, such as to capture the entire vortex flux. 

Additionally, eq. (|^) would allow for arbitrary world-surfaces of fluxes corresponding to even integer 
multiples of the elementary fluxes (|^). 
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to integer multiples of of course carry associated multiple field strengths, which would 
enter the topological density e^^i^^xF^fF^x- Such multiple field strengths were not allowed for 
in the requirement (|^) thus is not a new feature of the present investigation. The only 

freedom in (|^) is the choice of sign, i.e. the orientation of the vortex flux; this can be seen as 
a vestige of the color rotation freedom present in the characterization (|l|) which still remains 
after fixing the color direction of the gauge field via the Pauli matrix in (|^). In view of 



the nonorientability of generic vortex surfaces [T^, this residual freedom will be crucial 



This leads to a central issue, namely how the different individual cubes Ki are sewn 
together at the elementary squares (plaquettes) they share. By applying the specification (0) 
to a small loop encircling each plaquette on the lattice, a set of constraints on the choices of 
integers n{Ki), cf. eq. (j^), can be derived. To be precise, if the plaquette p extends from the 
lattice site x into the positive /x and u directions, and A, k denote the other two space-time 
directions]^, then consider the square loop C{p) given by the following sequence of corners, 

C{p) : y — ^ y + aex — ^ y + aex + ae^ — > y + ae^, — > y 

starting and ending at y = x + (a/2)(e^ + e^, — ca — e^), where a denotes the lattice spacing 
and Cp the unit vector in the p-direction. Clearly, then implies a constraint on the four 
cubes Ki whose boundaries contain the plaquette p (and which are thus intersected by C{p)), 
namely 

f dx^Ap = 

Jc{p) 

= [n{Kxix - ae^)) + n{K^{x)) - n{Kxix)) - n{K^{x - acx))] vra^ 

e {0,±7ra^} (8) 

where Kp{x) denotes the cube extending from x into all (positive) directions but the p 
direction. In physical terms, the plaquette p should carry at most one single unit of vortex 
flux (of either orientation). Every plaquette on the lattice generates one such constraint. It 
should be noted that this does not exclude any particular cube K on its own carrying an 
arbitrary integer n{K) in eq. (^; only from the point of view of choosing a gauge in which the 
gauge fields are as "smooth" as possible, cf. the general discussion in section ||, it is desirable 
to keep the magnitude of n{K) as small as can be achieved; how small this is in practice will 
be commented upon below. 

With the constraints (^ relating neighboring cubes Ki on the lattice, the global topology 
of the volumes S made up of such cubes enters the description. Namely, generic vortex 
surfaces S, and, concomitantly, volumes S spanning S, are not orientable . Formally, this 

""^""^A definite ordering of A, k is not necessary for tlic following construction. 
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expresses itself in the property that the constraints (Q) cannot all be simultaneously satisfied 
throughout space-time with a unique choice n{Ki) for all cubes Ki in the lattice. Rather, for 
non-orientable S, it is unavoidable to encounter frustrations in the attempt to satisfy all the 
constraints, i.e. one will find some plaquettes p within the volume S which yield the value 

dXf^A^ = ±27Ta^ , (9) 

C(p) 

thus carrying fiux double that of a vortex, cf. Fig. (left). Not being intended to appear as 
physical flux degrees of freedom within the random vortex surface model, these are the Dirac 
string world- surfaces which are unavoidable in a global, uniquely defined Abelian gauge field 
describing a non-oriented vortex surface; their boundaries correspond to monopole loops 



on the vortex sheets [^]. It is clear that such Dirac strings must arise when the vortex 
surface S is not orientable: In that case, there must be lines on the vortex surface where 
fiux orientation changes (the monopoles), and a globally defined Abelian gauge field, which 
satisfies continuity of fiux, must therefore contain Dirac strings supplying the change in vortex 
fiux at the monopoles. 

For reasons motivated in section 0, such (Abelian) Dirac strings are to be excluded from 
the gauge field in the present framework. Instead, the Wu-Yang construction will be used: 
Define the gauge field only locally on patches, sufficiently small such that all constraints 
(1) arising in any particular patch can be simultaneously satisfied. The different patches 
are then related via transition functions, cf. Fig. (right). The precise construction is as 
follows. The patches are the interiors of all the (2a) ^ hypercubes^ centered on the sites of 
the lattice, a denoting the lattice spacing. This is a natural choice, since these hypercubes 
will at the same time be the supports of the quark basis functions (finite elements) used in 
the construction of the Dirac operator matrix further below. In practice, these patches are 
sufficiently small to allow all the constraints (^ (and, of course, (^) arising on them to be 
satisfied. In order to achieve this, it furthermore turns out to be sufficient to allow for the 
values^ —2 < n{K) < 2 in (|^). An actual viable choice of n{Ki) for each of the 32 cubes Ki 
in a given patch is found by trial and error, starting from an initial cube, assigning values 



^^I.e., the cubes Ki making up the three-dimensional boundary of each hypercube in question are excluded 
from the patch. 

^■^As already hinted at above, this minimal choice of n(K) is motivated by the desire to formulate the gauge 
field in as "smooth" a gauge as possible, cf. the general remarks in section ^. On the other hand, it should 
not come as a surprise that, in general, one cannot be even more restrictive, i.e. the values n{K) = ±2 should 
indeed be allowed. Consider two parallel vortex surfaces of the same orientation near each other. Then the 
two segments of the volume E emanating from them may in general at some point (actually, a space-time 
surface) merge and run superimposed on one another from this merger surface onwards. Cubes Ki making up 
this coincidence volume must be associated with n{Ki) = ±2; the value n(Ki) = by contrast would imply 
a spurious double vortex (Dirac string) flux at the merger surface. In principle, the (gauge) freedom in the 
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Figure 1: Nonoriented gauge field support volume S in a two-dimensional slice of space-time. 
In such a slice, vortices show up as the solid points, volumes E emanating from them as the 
solid lines. The directions of the vortex fluxes at the boundaries of E, fixed a priori, here 
happen to be such that the segments of E emanating from the different vortices are forced 
to have mutually incompatible orientation. This is indicated by the arrows, which symbolize 
(directed) line integrals / A^rfx^; all of the line integrals indicated in the figures yield the 
same value, for definiteness -l-Trcr^. Left: Attempt at a global definition of the gauge field; 
the volume E then necessarily contains the indicated Dirac string, satisfying the property 
(^). Right: Wu-Yang construction, with space-time subdivided into two patches, on each of 
which the volume E is oriented. Gauge fields on the two patches are related by the transition 
function U . In the simple example depicted here, U can be chosen constant in the whole 
overlap region shared by the two patches. In generic configurations, the presence of other 
vortices nearby forces the space-time region of nontrivial U to be more localized, namely onto 
the immediate vicinity of E. This is discussed in more detail further below. 



to neighboring cubes such as to satisfy the corresponding constraints (H),®, and repeating 
until integers n{Ki) have been assigned to all cubes (note that the choice is in general not 
unique). In practice, this procedure uses an insignificant amount of computer time compared 
with the subsequent treatment of the Dirac matrix; thus, it is not necessary to invent more 
intelligent strategies for finding a viable set of integers n{Ki). 



spacc-tinic choice of S allows to deform S such as to avoid such a coincidence situation; in practice, when 
deformations are restricted by the coarse-grained, discrete lattice structure, this may not be directly possible. 
However, below, such n = ±2 volumes will indeed be reinterpreted as two n = ±1 volumes running parallel 
at a small distance from one another. In effect, this implies that, compared with its initial construction, the 
volume S as a whole may be augmented by additional (closed) volumes in order to consistently characterize 
the support of the gauge field A. 
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3.3 Matching patches - transition functions 



To facilitate a simple connection of the different patches via transition functions, two mod- 
ifications of the gauge field on each patch are still necessary. For one, as already hinted at 
above, any cube K which is associated with the value = 2 in at least one of the 

patches containing it is replaced by two cubes Ki and K2 parallel to K, and displaced from 
K hy a, small distance in the direction perpendicular to the cubes, cf. Fig. 0. The new cubes 
are assigned values n{Ki), n{K2) in all patches containing Ki and K2 as follows: 

n{K) = 2 ^ n{K{) = 1, n(K2) = 1 
n{K) = -2 ^ n{Ki) = -1, n{K2) = -1 (10) 
n{K) = ^ n{Ki) = 1, n{K2) = -1 

Clearly, the new gauge field on the two cubes Ki and K2 fulfils the constraint (0) defining 
vortex flux just like the original gauge field defined on K. All that has happened is that the 
phase picked up by any line integral intersecting K has been distributed onto the two cubes 
Ki and K2 which are slightly displaced from K. In effect, the splitting of K corresponds 
to a deformation of the volume S spanning the vortex world- surfaces S = (9S. The original 
E is augmented by closed volumes; such an operation does not change dT,, i.e. the physical 
content of the vortex configuration. Rather, it implies a {Z{2)) change of the gauge in which 
the gauge field is cast. As a result of the assignment ([T0|), only the values — 1 < n{K) < 1 
are still possible, at the expense of having introduced additional cubes K into the gauge field 
configurations. 



n=l 


n=l 




n=2 



n=l 



n=l 



n=l 



n=l 



Figure 2: Splitting of a cube K carrying n[K) = 2 into two parallel cubes Ki and K2 carrying 

The second change in the gauge field configurations which is still necessary stems from 
the following difficulty: In generic vortex configurations, one does not have the freedom to 
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choose constant transition functions throughout the overlaps between different patches, such 
as in Fig. |l| (right). Instead, the presence of other vortices nearby forces one to locahze 
the nontrivial parts of the transition functions to the immediate vicinities of the cubes K 
being rotated; in other words, the construction should be local to such an extent that the 
transition function can be specified independently for each cube K, cf. Fig. ^. However, such 
more sophisticated transition functions are space-time dependent and thus call for additional 
auxiliary gauge fields in some of the participating patches due to the inhomogeneous term in 
the gauge transformation law. Specifically, the modification consists in adding a pure gauge 
as follows. For any cube K associated with the value n{K) = — 1 on a given patch, consider a 
space-time region £{K) in that patch, bounded by two cubes and which are parallel 
to K, and which are displaced from by a small distance e, as well as a small volume Kq (of 
extension e in one of its directions) connecting Kj^ and i^-, such as depicted in Fig. ^. Note 
that any definition of £{K) such that £{K) contains K, but no points more distant from K 
than e, serves the purpose of the present construction. If desired, it is quite admissible to 
deform Kq in a way which excludes from £{K) some points which are closer to K than e. 
Note also that the distance e is intended to be even smaller than the distance between two 
cubes Ki and K2 used in the cube-splitting procedure discussed in the previous paragraph. 
Consider furthermore the gauge transformation 

-ta^ for xe8{K) ^^^^ 

and the gauge field a{K) = iUlf^j^-^dUsi^x) induced by U£(^k), which has support on the bound- 
ary of S{K) described above. 

On any patch on which n{K) = — 1, the field a is added to the gauge field in that patch, 

A(K) I = 1 (12) 

\A{K) + a{K) for n{K) = -1 ^ > 

for each cube K in all patches. This modification does not change any Wilson loops nor any 
line integrals (0) taken in the configuration; also the physically relevant field strength content 
of the gauge field on the patch in question is not changed, apart from the following subtlety: 
In general, it is not excluded that Kq may be forced to intersect the volume S which spans 
the vortex world-surfaces S = (9E, and which supports a nontrivial gauge field A, cf. Fig. |[ 
In this case, a spurious flux is generated on the intersection surface through the commutator 
term in the field strength, Ff^i, oc [A^,ay]. Thus, in some cases, spurious fluxes of the Dirac 
string type are reintroduced; however, the crucial point is that these fluxes have now been 
cast in a form in which infrared quark propagation manifestly is not influenced by them. 
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Figure 3: Space-time patch with three cubes K, Ka and Kb which are part of the gauge 
field support volume S; in the two-dimensional slice displayed, these cubes show up as the 
solid lines. For definiteness, n{KA) = n{KB) = 1 and n{K) = —1; there is a vortex at the 
center of the patch. On a neighboring patch which also contains K, this cube shall carry the 
value n'{K) = 1. The region on which the associated transition function U is nontrivial is 
localized as specified in the figure; this simultaneously implies supplementing the gauge field 
in the depicted patch with an auxiliary field d = iU^dU, see text. The field d has support 
on the volumes K^, K_ and Kq, which is where U changes as a function of space-time. In 
general, it cannot be avoided that Kq intersects parts of S (in the absence of Ka and Kb, one 
could still deform Ko such as to circumscribe K). Such intersections imply spurious fiuxes, 
see text. 



Namely, consider any coupling matrix element 

J d'^xql dq2 (13) 

taken between smooth]^ quark wave functions qi and g2- Such an integral is negligible, since 

^^The alert reader may notice that, in general, one will also have to consider the case where one of the quark 
wave functions is discontinuous due to the presence of a nontrivial transition function e.g. q2 = Uq'2 

with smooth Here, represents a smooth quark field on an adjacent space-time patch, which must be 
gauge-transported to the current patch via the transition function U for the purpose of evaluating the matrix 



element. This discussion is deferred to the beginning of section 4.1, where the quark basis functions are 
defined, in order to be able to treat the issue in a more definite manner. The statement that integrals of the 
form ([13) are negligible will indeed remain true even in this more general case. 
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the fields on and differ precisely by a sign, and and are located only a very 
small distance e apart; furthermore, the contribution from Kq is of order e due to the small 
extension of Kq in one of its space-time directions. Thus, the smooth quark modes of an 
infrared effective theory are insensitive to the supplementary field a. It is not excluded that 
a more sophisticated construction exists which avoids the residual spurious field strengths 
induced by the additional presence of d altogether; however, in view of the insensitivity of 
infrared quark propagation with respect to a, it does not appear necessary to develop such 
a construction at this stage. In a sense, the situation is opposite to the case of the Dirac 
string fluxes which one would obtain by naively defining a global Abelian gauge field such 
as discussed further above; in that case, the fact that one is working with a truncated quark 
basis precludes the solutions of the Dirac equation adjusting such as to make the Dirac strings 
invisible. In the present case, it is precisely the infrared sector of smooth quark wave functions 
which manifestly decouples from any residual spurious fluxes introduced by the addition of 
the field d. Note, however, that this is one of the main points preventing a straightforward 
generalization of the present gauge field definition to explicitly thickened vortex fluxes and 
gauge field volumes. Such a generalization would entail finite e, and as a consequence, the 
spurious additional field strengths induced by adding the field d as discussed here could no 
longer be argued to decouple. Thus, a more sophisticated construction of the gauge fields on 
each patch and the associated transition functions would need to be specified. 

Note furthermore that, at this point, one still has a certain amount of freedom as regards 
the global definition of the gauge field, which will be discussed in more detail below in 



section p.4| . In particular, it is still possible to vary the space-time density of nontrivial 
matchings between neighboring patches, i.e. nontrivial transition functions. Comparing the 
results obtained with constructions which mutually differ in this respect will give an indication 
of the influence of the rough disorder in the gauge field associated with the presence of those 
nontrivial transition functions. 

Having prepared the gauge field configurations on each space-time patch as described 
above, it is straightforward to give the corresponding transition functions between patches. 
Consider two patches, distinguished in the following by denoting all quantities in one patch 
with a prime and their counterparts in the other patch without a prime. Consider furthermore 
the space-time region in which the two patches overlap and the set of cubes Ki located 
there, associated with values n{Ki) and n'{Ki) in the two patches. The associated transition 
function in that overlap region then is 

f/=n^i^S^""'^'''^'^' (14) 

with Us(^K) given by (0), i.e. the gauge fields A (including, if present, the auxiliary field d) 
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and the quark fields q in the two patches are related by 

A' = UAU^-iUdU^ (15) 
q' = Uq (16) 

3.4 Options in the gauge field construction 

The gauge field construction presented above still allows for different model options in two 
respects. In the numerical work further below, a variety of these options will be explored 
in order to be able to appreciate the dependence of the results on the gauge field modeling, 
which, after all, contains a range of ambiguities and idealizations as presented in the preceding 
sections. 

One of the options one still has lies in the choice of the volume S spanning the vortex 
world-surfaces. A particular such volume S, the construction of which was repeatedly deferred 
above, can be obtained as follows. As discussed in detail in |^, Monte Carlo updates of 
vortex configurations in the present model are always carried out simultaneously on all six 
plaquettes making up the surface of an elementary three-dimensional cube in the lattice; this 
serves to keep the vortex surfaces closed as they are being updated. As a result, a volume 
S whose boundary defines the vortex surfaces can be trivially obtained by keeping track of 
the cubes being updated. To be precise, if one starts with an empty lattice, (f){Ki) = 1 for 
all three-dimensional cubes Ki in the lattice; whenever an update of (the surface of) a cube 
K is subsequently accepted, one concomitantly changes (p{K) —(j){K). In this way, at any 
point in the calculation, a viable E, made up of the set of with (f){Ki) = —1, is known. 

On the other hand, this is not the only possible S. Instead of working directly with 
the above S, which will generally be rather random and rough, it is possible to first apply a 
smoothing algorithm to S. Since the only physically relevant information lies in the boundary 
of E, it is admissible to add to E, or take away from E, arbitrary closed three-volumes. In 
practice, this implies sweeping through the lattice, considering each elementary four-cube 
in turn, and simultaneously updating all eight three-cubes making up its boundary, i.e. 
4>{Ki) —(f){Ki), whenever this leads to a reduction in the total volume of E. The two 
options for E are related by a {Z{2)) gauge transformation which however does not vary 
smoothly in space-time. Consequently, invariance of the Dirac spectrum is not guaranteed 
in a calculation using a truncated quark basis, as already mentioned in section Thus, 
the two options represent different possible models which will be explored in the numerical 
work below; the options will be referred to as "random E" and "smooth E", respectively. 
Of course, under the aspect of formulating the gauge field in as smooth a gauge as possible 
in order to minimize artefacts stemming from rapidly varying gauge transformations, the 
smooth option seems the preferred one. 
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The second option one still has in the construction of the gauge field is the following: 
Above, the gauge field A is initially constructed independently on each space-time patch. 
As a result, the relative color orientation of adjacent patches, ultimately encoded in the 
transition functions on the respective overlap regions, is random; there is a high density of 
overlap regions supporting nontrivial transition functions. However, for any given patch, the 
transformation A — > —A leads to an equally admissible gauge field on that patch; this simply 
corresponds to a gauge transformation U = —ia"^ constant throughout the patch. Thus, 
after the initial construction of the gauge field A (but before the addition of the auxiliary 
gauge field a discussed in section one has the alternative of performing sweeps through 
the lattice in which one considers a transformation of the gauge field A —A on each 
individual patch in turn. Acceptance of such a transformation can be biased such that as 
many as possible of the resulting transition functions to neighboring patches become trivial. 
In other words, the transition functions are made maximally smooth. In the language of 
globally defined gauge fields with monopole and Dirac string singularities, this corresponds 
to minimizing the space-time density of such singularities. Accordingly, the option of fixing 
the gauge field orientation at random on the patches will below be referred to as the "ran- 
dom monopole" gauge, whereas the option of maximally smooth transition functions will be 
referred to as the "minimal monopole" gauge. Note that such an alternative, either choosing 
the monopoles at random or minimizing their density, was already considered in 0,[14|,[26 . 
The topological susceptibility turned out to be virtually independent of the monopole den- 
sity, a fact which could be understood in terms of the space-time properties of generic vortex 



surfaces 0,[14|,|26l. 



4 Vortex- Quark coupling 
4.1 Quark basis and Dirac matrix 

As already indicated further above, the Dirac equation will be solved using the finite element 
method. This means that the truncated quark basis is constructed using space-time wave 
functions f{x — xq) associated with, and localized around, each lattice site xq (the Dirac and 
color structure will be supplemented further below). The wave function / is a product over 
the four space-time directions, 

f{z) = f[h{z,) (17) 

i=l 
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and h is piecewise linear, 



{t + a for -a < t < 
a - t for < t < a (18) 
else 

where a denotes the lattice spacing. To be precise, the space-time wave function associated 
with the lattice site xq takes the form f{x—xo) on the patch which is simultaneously associated 
with xq. Viewed from other patches, this form is of course modified by the corresponding 
transition functions^. An exception to the above form occurs in the case when xo is located 
on the boundary where the space-time torus is sewn together, i.e. when at least one of the 
vector components of xq, say xoj, vanishes. In this case, the corresponding factor(s) h{zj) in 
(0) is (are) replaced by 

{—a — Zj for —a < zj < 
a — Zj for < Zj < a (19) 
else 

In other words, the quark wave functions are forced to obey antiperiodic boundary conditions 
in all directions. Physically, this is only really necessary in the Euclidean time direction; the 
space directions on the other hand should always be taken sufficiently large such that bulk 
physics does not depend on the behavior at the boundaries. Note that, as far as the torus 
periodicity is concerned, the gauge fields discussed further above obey periodic boundary 

^^Note that this means that quark basis functions are only strictly smooth on the patch they are defined 
on, whereas, viewed from the color frame associated with a neighboring patch, they in general exhibit non- 
smooth behavior near selected three-dimensional cubes K; after all, the corresponding transition functions 
can vary rapidly on the very small length scale e entering the definition of £{K) in As one moves on 

such a cube K from one lattice site to the next, a general quark wave function, linearly combined from the 
basis functions, continuously interpolates between a form which is strictly smooth on the patch associated 
with the initial site to a form which is strictly smooth on the patch associated with the final site. It should 
be emphasized that even with this generalized (not strictly smooth) ansatz for the quark wave functions, the 
statement made in connection with eq. (p^ ) remains true, i.e. the auxiliary field a cancels in all Dirac matrix 
elements and in this sense does not influence quark propagation. This follows from the property that the 
part of the wave function rapidly varying around a given cube K, viewed as a function of the coordinate 
perpendicular to iiT, is symmetric around K. As a consistency check, if one considers the free Dirac equation 
in patched space-time with nontrivial transition functions and the corresponding auxiliary fields a, the Dirac 
matrix is manifestly identical to the one with trivial transition functions. In the framework of the Wu-Yang 
construction, the transition functions are nothing but a part of the choice of gauge; the fact that they (and, 
as a consequence, the quark wave functions) are not strictly smooth is merely a specific instance of the model 
compromise already discussed in section ^ namely that the fields in the present construction are not given 
in the smoothest possible gauge, but that the gauge is partly dictated by considerations of manageability. 
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conditions; translation along one of the torus directions by its length leads back to the same 
patch, with the same gauge field. 

A further comment is in order concerning space-times of extent a in the Euclidean time 
direction. While the vortex surfaces and the associated gauge fields can be straightforwardly 
defined on such a space-time, the truncated quark basis should include the two lowest Mat- 
subara frequencies, which are degenerate, on an equal footing. Thus, in practice, while the 
vortex configurations are still generated on the original lattice with spacing a, these configu- 
rations are subsequently copied onto a lattice with two spacings a/2 in the time direction in 
order to allow for an (approximate) accomodation of both of the basis functions corresponding 
to the lowest Matsubara frequencies. 

The above set of wave functions / is supplied for each color and Dirac component, i.e. 
the complete quark basis can be labeled as 

g(6, i, xo) = Cbdif{x - xq) (20) 

where q, is the 6th unit vector in (two-dimensional) color space, di is the zth unit vector 
in (four-dimensional) Dirac spinor space, and xo is a site on the space-time lattice. Note 
that quark wave functions are linearly combined from these basis vectors using, in general, 
complex coefficients. In this basis, it is straightforward to calculate analytically the matrix 
elements of the (Euclidean) Dirac operator 

P = - A^) (21) 

using the gauge field A constructed in section ^. This becomes particularly simple when 
the diverse small lengths introduced in section |3.3| are in fact taken to be negligibly small, 
specifically, the distance separating the two parallel cubes Ki and K2 in eq. (0), and the 
(even smaller) width e of the domain S{K) used in eq. (|ll|). The Dirac matrices were in 
practice taken in the chiral representation. 
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(22) 

When evaluating matrix elements between states defined on adjacent patches, of course the 
corresponding transition function must be used to transport all functions involved to a definite 
color frame, i.e. for unequal Xo,Xq, 

{b'i'xQ\p\bixo) = Cb'di' J d'^x f{x - XQ)U{x;xQ,Xo)pf{x - XQ)cbdi (23) 

where U{x;Xq,Xo) denotes the transition function connecting the patches centered at xq and 
x'q, the assignment of primed and unprimed patches having been chosen as in eqs. ([T^)-([16|). 
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Note that the quark basis used here is not orthogonal, i.e. the overlap matrix M made 
up of the matrix elements 

{h'i'x'^\UxQ) (24) 

is not the unit matrix. Thus, to obtain the eigenvalues A„ of the Dirac operator, one must 
solve the matrix equation 

pZn = XnMZn (25) 

Furthermore, since the square of the Dirac operator and the overlap matrix M are block- 
diagonal in the Dirac indices, in practice it is advantageous to instead diagonalize the squared 
Dirac operator, i.e. to consider the equation 

M-^D^M-^Dzn = \\zn (26) 

which is only two-dimensional in Dirac spinor space, D denoting the matrix of the operator 

b = idQ- Aq- ia\idi - Ai) (27) 

Note that the Dirac matrix constructed here is manifestly Hermitean and chirally symmetric, 
i.e. its eigenvalues are real, and when they are nonzero, they trivially appear in pairs of 
opposite sign. On the other hand, since the quark basis is finite-dimensional, no exact zero- 
mode solutions occur. As a consequence, the spectrum obtained here by construction strictly 
only contains the non-zero mode sector of vanishing chirality. This is different from lattice 
Dirac operators, in which both chiral symmetry is broken and the quark basis is truncated; 
the continuum limit then corresponds to a combined restoration of chiral symmetry and basis 
completeness - this in principle also allows one to capture exact zero modes in the limit[^ 
By contrast, in the present continuum approach, the limit of exact chiral symmetry is taken 
already at fixed basis truncation. With this order of limits, chiral modes cannot be recovered 
and must be obtained separately by different means. Such an investigation, e.g. via the 
spectral flow of the Dirac operator in the presence of an explicit chiral symmetry-breaking 
parameter, lies outside the scope of the present work. For the purpose of describing the 
chiral condensate, the exact zero mode sector is irrelevant in the infinite volume limit, simply 
because the topological susceptibility {Q"^) /V is a finite quantity {Q denotes the topological 
charge, each unit of which is expected to be associated with a quark zero mode via the 
Atiyah-Singer index theorem, and V is the space-time volume under consideration). As a 
consequence, the contribution of the exact zero modes to the total Dirac spectral density 
becomes negligible as \^ ^ oo. 

Ultimately, this is connected with the fact that the standard lattice cutoff regularizes the theory gauge- 
invariantly in the ultraviolet, as opposed to the fixed momentum cutoff implied by an infrared effective 
framework such as discussed here. 
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4.2 Remarks on topology and torus twist 

Before applying the construction presented in the previous sections to the configurations of 
the random vortex surface model, a comment on their topology is in order. The topological 
charge of these configurations is quantized in half-integer units 0. Usually, such half- integer 
charge is thought to require twisted boundary conditions for the gauge fields on the torus 
These are boundary conditions of the type 

+ = (28) 

where Li is a vector in the i-direction of the length of the torus in that direction, and Ui 
denotes a gauge transformation. The latter must obey the consistency condition 

Ui{x)Uj{x + L,) = ±U,{x)Ui{x + L,) (29) 

in order to give a unique, well-defined relation between A{x) and A{x + Li + Lj). The case of a 
minus sign on the right hand side of ( p9D is referred to as a twisted boundary condition^^ In 
particular, fundamental quark fields cannot be defined on a torus with such twisted boundary 
conditions; applying the gauge transformations on either side of ( pOj) to a quark field q{x) 
implies 

q{x + Li + Lj) = -q{x + L, + Lj) , (30) 

i.e. q = 0. 

In contrast to this, the construction presented in the previous sections uses purely pe- 
riodic boundary conditions for the gauge fields on the torus, and in particular also allows 
the definition of (antiperiodic) quark fields. With the choice (|T^) of transition functions. 



the quark fields are uniquely defined on every patch and no inconsistency such as (^) can 
arise. The construction includes arbitrary vortex world-surface configurations generated in 
the random vortex surface model, in particular ones with half-integer topological charge. 

This does not contradict the usual statement that half-integer topological charge on a 
torus requires twisted boundary conditions, because the space-time manifold used in the 
construction is not a torus in the strict sense. Rather, the spurious fiux surfaces induced 



^^In the vortex language, the mechanism by which twisted boundary conditions allow for half-integer 
topological charge takes the following form. Twisted boundary conditions permit the existence of a single 
vortex world-surface corresponding to a two-dimensional plane in the torus. On its own, such a plane does 
not represent the boundary of any three-dimensional volume, and thus cannot occur for standard periodic 
boundary conditions. In the presence of two such planes which are completely perpendicular to each other, 
there is precisely one intersection point of these planes on the torus, which represents the lone contribution 
to the topological charge ||],[^, of modulus 1/2. Note that the vortex surface configuration in this case 
is oriented; by contrast, for compact vortex surfaces which do not rely on specific space-time boundary 
conditions, nonvanishing topological charge requires non-orientedness |23]. 
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by the auxiliary fields a discussed in section p.3| , which in general cannot be avoided, must 
be excised from space-time, leading to a complicated, multiply connected structure[^. By 
contrast, the case usually referred to as a twisted torus corresponds to the specific case where 
there are no such excisions, and there is only a single patch encompassing all of space-time, 
merely with nontrivial transition functions, cf. eq. (0), where opposite ends of the patch are 
sewn together to make up the toroidal space-time structure. This represents a very rigid and 
direct relation between the global shape of space-time and the patching, and the combina- 
tion of the two leads to the constraint (pQ]). On the other hand, in the present more flexible 
construction, the patching is completely divorced from the global torus structure; already the 
local structure of space-time, on scales related to the strong interaction scale, is topologically 
nontrivial in a dynamically determined fashion. In this setting, the additional toroidal peri- 
odicity constraints are not as crucial in determining the global topological properties of the 
vortex configurations as in the one-patch, twisted torus case (nevertheless, these periodicity 
constraints are still there, with the associated finite size effects). 

A particular example of a vortex configuration illustrating this point is given in Fig. ^ 
This configuration has topological charge (5 = 1/2 (for details on evaluating the topological 
charge of a vortex configuration, cf. 0). On the other hand, this world-surface can be 
represented as the boundary of a three-dimensional volume which is localized in the region 
of the vortex. Thus, within the present construction, this configuration can be encoded in a 
gauge field A with support in a compact domain, and A = in all directions if one moves 
away from this domain. Therefore, this configuration can be defined on space-times of any 
shape, without relying on special assumptions concerning the behavior of the gauge fields 
at the space-time boundaries. However, the surface is non-orientable; a possible position of 
the associated (non-contractible) monopole loop on the vortex world-surface is given by the 
thick line in Fig. ^. This implies that, within the present construction, a surface spanning the 
monopole loop, carrying spurious nonvanishing field strength, must be excised from space- 
time. Concomitantly, the configuration cannot be defined on a single space-time patch; it 
requires at least two space-time patches, related by nontrivial transition functions. In this 
way, the configuration induces a nontrivial space-time topology within the region occupied 
by it; the "twist" permitting topological charge 1/2 is localized in that region, instead of 
being encoded in nontrivial (twisted) torus boundary conditions. 

At this point, there seems to be no apparent reason to exclude the possibility of an alter- 
nate gauge field description of such configurations which indeed transfers the twist required 
for half-integer topological charge into twisted torus boundary conditions. Note e.g. that the 

^^Note also that, in such a multiply connected, patched space-time, it is a priori not clear that the connection 
between topological charge and zero modes of the Dirac operator, i.e. the index theorem, takes the same simple 
form as in a standard one-patch space-time. 
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nonabelian theory allows for pure gauges which, in their diagonal color components, contain 
Dirac monopole loops spanning open Dirac string world- surfaces (by contrast, in the Abelian 
gauge theory, only closed Dirac strings are pure gauges). An explicit example of such a 



pure gauge is given in p3|; it contains long-range nonabelian gauge fields. Thus, one could 
envisage using such gauges to cancel any spurious Dirac string type fluxes in a gauge field 
description of the configuration depicted in Fig. ^, at the expense of introducing nonvanishing 
gauge fields at large distances; these in general will be influenced by the space-time boundary 
conditions. 




Figure 4: Sample closed vortex surface configuration. At each lattice time wq, shaded plaque- 
ttes are part of the vortex surface. These plaquettes are furthermore connected to plaquettes 
running in time direction; their location can be inferred most easily by keeping in mind that 
each link of the configuration is connected to exactly two plaquettes (i.e. the surface is closed 
and contains no intersection lines). Note that the two non-shaded plaquettes at uq = 2 are 
not part of the vortex; only the sets of links bounding them are. These are slices at no = 2 
of surface segments running in time direction from no = 1 through to no = 3. Sliced at 
no = 2, these surface segments show up as lines. The surface is non-orientable; the thick line 
at no = 2 indicates a minimal non-contractible (monopole) loop on the surface at which the 
orientation is inverted. Note also that the topological charge Q = 1/2 of the configuration is 
carried purely by surface writhings; there are no surface self-intersection points. The author 
gratefully acknowledges R. Bertie and M. Faber for providing the MATHEMATICA routine 
with which the image was generated. 



Note furthermore that even the fact that a twisted torus does not directly permit a 
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definition of fundamental quark fields should be seen more as a practical disadvantage than 
a fundamental one. By piecing copies of a twisted torus together such as to make up a 
larger torus with appropriate integer multiple extensions in the four space-time directions, 
one arrives at a description in terms of gauge fields with periodic boundary conditions on 



the larger torus. On this torus, quark fields can then be straightforwardly defined p8 



The price one pays consists in having to work on a larger space-time manifold; of course, 
from the practical point of view of numerical Monte-Carlo calculations, this can be a decisive 
disadvantage. 



5 Quenched Dirac spectrum of the random vortex sur- 
face model 

In the following, results for the quenched ensemble average of the Dirac spectrum in the 
framework of the random vortex model of infrared Yang-Mills theory will be presented. 
The dynamics determining this model ensemble and its physical interpretation are discussed 
at length in [Q. The model describes closed vortex world-surfaces composed of elementary 
squares (plaquettes) on a hypercubic lattice. An ensemble of such world-surfaces is generated 
using Monte Carlo methods, where, in order to keep the surfaces closed, an elementary update 
of a vortex world-surface configuration always affects all six plaquettes making up the surface 
of a three-dimensional cube in the four-dimensional lattice; any of those six plaquettes which 
were not part of a vortex surface before the update become part of a vortex surface after 
the update, and vice versa. As mentioned in section p.4| , keeping track of the updated three- 
cubes directly allows the construction of the three-volume S which is needed as input for 
the definition of a corresponding vortex gauge field, cf. section |3.1| . The acceptance criterion, 
i.e. the Boltzmann factor exp(— S"), for any given update depends on the curvature of the 
vortex surfaces. An additive action increment c = 0.24 is contributed to S by each pair of 
vortex plaquettes on the lattice sharing a link, but not lying in the same plane; the value 
of c was chosen such as to quantitatively reproduce the confinement properties of SU (2) 
Yang- Mills theory [|l| . This vortex ensemble simultaneously yields a result for the topological 
susceptibility which quantitatively agrees with full SU{2) Yang-Mills lattice calculations 

Given a particular vortex configuration, via the associated three- volume E, the Dirac 
matrix (in practice, a 2 x 2 Dirac block of its square, cf. eq. (^)) can be constructed as 
described in sections ^ and |4.1| . The 30 lowest-lying eigenvalues of this matrix were obtained 
using the ARPACK package ||3T| , and (their positive square roots) binned according to their 



magnitude, yielding the density of eigenvalues (per unit space-time volume) p(A) of the Dirac 
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operatoi]^. The behavior of this distribution near A = is related to the chiral condensate 
{qq) via the Casher-Banks formula |]32[|,p3 

{qq) = - Jim 2m H d\ -7rp(0) (31) 

where m denotes the quark mass and V the space-time volume. A range of results obtained 
in this way is displayed in Figs. |^ and ^, with emphasis on the variation with volume. The 
two main features of these eigenvalue distributions are a bulk behavior which extrapolates to 
a finite non-zero value at A = 0, and an additional enhancement for very small eigenvalues 



which presumably signals the onset of divergent behavior at A = 0. In view of (|3T|), such a 
divergence implies a diverging chiral condensate in the chiral limit m 0; at finite quark 
mass m, the enhancement near A = generates an increase in the chiral condensate compared 
with the value one would extract from the bulk of the spectrum. On the other hand, in (|3T|) 
the infinite volume limit is to be taken before the chiral limit, and Figs. |^ and ^ suggest that 
the enhancement of the spectra near A = weakens as the space-time volume is increased. 
However, the limited statistics of the measurements precludes a solid extrapolation of the 
aforementioned enhancement to infinite volume; it is unclear whether it completely disappears 
as — oo or whether it approaches a nontrivial limit. 

Divergences in the quenched chiral condensate reminiscent of the behavior observed here 
have been noted in a variety of contexts, and both of the aforementioned alternatives concern- 
ing the ^ oo limit have been discussed. An ensemble of vortex configurations extracted 
from the full SU{2) lattice Yang-Mills ensemble using an appropriate gauge fixing and projec- 
tion procedure yields a divergent chiral condensate in the chiral limit This phenomenon 
is also observed in quenched lattice calculations using domain wall fermionsQ [0 . The main 
effect in this case is due to exact Dirac zero modes resulting from a nontrivial topological 
charge of the gauge configurations; this effect vanishes for — > oo as has been argued at the 



end of section [4.1| . On the other hand, this effect may also hide more subtle divergences which 
persist in the infinite volume limit. Such more persistent divergences have been argued to 
appear both in quenched chiral perturbation theory and in instanton models In the 



latter models, the strength of the divergence strongly decreases as one leaves the dilute in- 
stanton gas limit and allows the associated quark zero modes to interact strongly. Note that, 
in contrast to the instanton picture, the present vortex picture does not possess an obvious 
dilute limit with weakly interacting quark quasi-zero modes. In view of this, the observation 



^^Strictly speaking, this yields the part of the spectral distribution located on the positive eigenvalue axis. 
However, due to manifest chiral symmetry, the distribution as a whole is symmetric about eigenvalue zero; 
therefore, the distribution on the negative eigenvalue axis does not need to be considered explicitly. 

^°Using more conventional lattice Dirac operators, such behavior is not detected, presumably being veiled 
by lattice artefacts pO|. 
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of a persistent divergence in the Dirac spectrum as — > oo in the instanton model does 
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Figure 5: Dependence of Dirac spectral density p on the extension of space-time. Left: Gauge 
field constructed using random monopoles and random S; data symbolized by x corresponds 
to a 3^ lattice, crosses correspond to a 4^ lattice, and circles to a 5^ lattice. The reader is 
reminded that the lattice spacing is 0.39 fm. Right: Gauge field constructed using minimal 
monopoles and random E; data symbolized by x corresponds to a 3^ x 2 lattice, crosses 
correspond to a 4^ x 2 lattice, and circles to a 6^ x 2 lattice. Note that an extension of 
the lattice in Euclidean time direction of two spacings is equivalent to a temperature of 
T = 0.83 Tc. 

not straightforwardly imply the same qualitative effect in the vortex picture. More detailed 
information on this question may result from a consideration of the space-time structure 
of low-lying quark modes in the vortex model; such an investigation lies beyond the scope 
of the present work. Note also that, in unquenched calculations, the spectral density near 
zero eigenvalue will be dynamically suppressed due to the weighting with the determinant of 
the Dirac operator. Nevertheless, the aforementioned instanton model calculations still yield 
divergent behavior at A = 0, which however disappears as m — > such as to yield a finite 
chiral condensate in the chiral limit. 
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Figure 6: Dependence of Dirac spectral density p on the extension of space-time. The 
gauge field is constructed using minimal monopoles and smooth E. Data symbolized by x 
corresponds to a 3^ x 2 lattice, crosses correspond to a 4^ x 2 lattice, circles to a 5^ x 2 lattice, 
open boxes to a 6'^ x 2 lattice, and filled boxes to a 7^ x 2 lattice. 



If the Dirac spectrum retains divergent behavior at A = even in the infinite volume 
limit, then the issue of the phenomenological relevance of this behavior becomes a subtle 
quantitative question. In view of (|HTD, for sufficiently large quark masses the extra contri- 
bution to the chiral condensate resulting from the divergence at A = becomes negligible 
and the condensate is dominated by the bulk spectral density. The question then remains 
whether realistic quark masses are large or small in this sense. While this caveat must be 
kept in mind, further below, numerical values for the chiral condensate will be quoted which 
correspond to the A = extrapolation of the bulk of the spectral density, disregarding the 
enhancement of the spectrum very near A = 0. 

Turning now to the properties of the bulk of the spectrum in more detail. Fig. ^ displays 
results for the "random S" construction of the gauge field with either minimal or random 
monopoles, respectively. In these cases, the bulk value evidently has already converged well 
to the large volume limit at a spatial extension of the universe corresponding to 4 lattice 
spacings, i.e. 4a = 1.56 fm. In Fig. |^, exhibiting results for the "smooth E" construction 
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with minimal monopoles, the convergence is not quite as fast; there is still a discrepancy of 
the order of 10 % between the universe with spatial extension 1.56 fm and the universe with 
spatial extension 1.95 fm, the results in the latter case being indistinguishable from those 
on larger lattices^. The finite size effects were also considered in a number of other cases: 
Random S with random monopoles on x 1 lattices^, random E with minimal monopoles 
on lattices and on x 1 lattices, and smooth S with minimal monopoles on lattices. 
All these cases displayed a convergence to the large volume limit comparable to or better 
than the ones exhibited here explicitly. Also in the smooth S case with random monopoles, 
the bulk of the spectrum on a 3^ lattice turns out to be indistinguishable from the one on a 
4^ lattice. 

In the following, Dirac spectra will be shown for different choices of the gauge field con- 
struction and different temperatures at a spatial extension of the lattice of 4a = 1.56 fm. 
While one should be aware that this extension in some cases still allows for finite size effects 
of the order of 10 %, it does not suffer from the aforementioned subtlety (cf. footnote) arising 
for lattices with odd numbers of lattice sites in at least one of the space-time directions, and 
it also allows for superior statistics compared with extensions corresponding to the next even 
multiple of the lattice spacing, 6a = 2.34 fm. The results are depicted in Figs. 0-^, and also 
summarized in Table |l]. 

Evidently, in the confined phase, all models yield a bulk spectral density which extrapo- 
lates to a nonzero value at A = 0, thus signaling the spontaneous breaking of chiral symmetry. 
Quantitatively^, the zero-temperature quenched chiral condensate obtained for the different 
gauge field constructions varies between -0.1 GeV^ and -0.175 GeVl For comparison, the 
value obtained in lattice Yang-Mills calculations ||20| corresponds to —0.09 GeV^ (with the 

^^Note that the approach to the infinite volume limit is not monotonous. This is presumably connected 
to the following subtlety: In the case of a lattice with an odd number of lattice sites in any one of the 
space-time dimensions, the space of quark wave functions used here contains a mode which is not annihilated 
by the derivative operator, but whose image after acting with the derivative operator has no more overlap 
with any quark wave function. Consequently, the free Dirac matrix then contains corresponding spurious 
zero eigenvalues; as a result, the free Dirac equation is only solved correctly if one uses an even number 
of lattice sites in all directions. In the presence of an interaction with the gauge field, the aforementioned 
problem of vanishing overlap largely disappears, as is evidenced by the similarity of the spectra obtained for 
the cases of either all even or some odd numbers of lattice sites in the four space-time directions; however, 
the non-monotonous approach to the infinite volume limit is presumably a vestige of the subtle difference 
between these two cases. 

^^The reader is reminded that, in the special case of a space-time of extension a in the time direction, a 
subdivision of the lattice in that direction is introduced when defining the quark basis, in order to accomodate 



both of the lowest Matsubara modes, cf. section 4.1 



^•^The scale in the random vortex surface model is fixed by equating the zero-temperature string tension 
with (440 MeV)^, simultaneously implying a lattice spacing of a = 0.39 fm in the model with curvature 
coefficient c = 0.24, cf. |]. 
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Figure 7: Dependence of Dirac spectral density p on temperature T, at a spatial extension 
of the lattice of 4a = 1.56 fm. Left: Gauge field constructed using random monopoles and 
random S; data symbolized by x corresponds to T = 0, circles correspond to T = 0.83 Tc, 
and crosses to T = 1.67Tc. Right: Gauge field constructed using minimal monopoles and 
random S; data symbolized by x corresponds to T = 0, circles correspond to T = 0.83 T^, 
and crosses to T = 1.67 T^. 



zero-temperature string tension fixed to the same value as in the present model). In particu- 
lar, both of the random S models come very close to the lattice Yang-Mills result. However, 
it should be kept in mind that the actual value of {qq) does not have direct phenomenological 
meaning, since it is not a renormalization group invariant quantity^. Only after multiplying 
with the quark mass, one obtains a physically meaningful quantity which enters e.g. current 
algebra relations. In this sense, a measurement of {qq) can be viewed as representing no more 
than an accessory to defining the current quark masses within the model. However, the fact 

^''At this stage, the relation between the ultraviolet cutoff scheme of the present vortex model on the one 
hand and, on the other hand, standard renormalization schemes such as the ones used on the lattice or in 
perturbation theory is not fixed. In principle, matching different schemes can lead to large correction terms; 
this is well known in the case of matching lattice renormalization with standard perturbation theory schemes, 
where e.g. the chiral condensate acquires corrections of the order of 30% pO|. 
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Figure 8: Dependence of Dirac spectral density p on temperature T, at a spatial extension 
of the lattice of 4a = 1.56 fm. Left: Gauge field constructed using random monopoles and 
smooth E; data symbolized by x corresponds to T = 0, circles correspond to T = 0.83 Tc, 

and crosses to T = 1.67 T^. Right: Gauge field constructed using minimal monopoles and 
smooth E; data symbolized by x corresponds to T = 0, circles correspond to T = 0.83 T^, 
and crosses to T = 1.67 Tc. 



that the vortex model yields values for the chiral condensate which are not separated from 
lattice values by unnaturally large factors gives rise to the expectation that this model suc- 
ceeds in capturing the relevant physics leading to spontaneous chiral symmetry breaking in 
the confined phase of the strong interaction. More stringent quantitative tests of low-energy 
quark physics within the vortex model demand the calculation of additional observables; this 
is deferred to future work. 

Turning now to the behavior of the chiral condensate as one crosses the finite-temperature 
dcconfinement transition, the different models evidently display much larger, qualitative, 
differences. The random S, random monopole model yields a considerably enhanced chiral 
condensate in the deconfined phase; the other models show a reduction of the condensate 
above the critical temperature. However, this reduction is rather slight, except in the case of 
the smooth E, minimal monopole model. This latter case is the only one which qualitatively 
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Figure 9: Detailed plot of the (third root of the) Dirac spectral density p obtained at a 
temperature of T = 1.67 T^. for a spatial extension of the lattice of 4a = 1.56 fm using the 
minimal monopole, smooth S construction of the gauge field. The chiral condensate does 
not vanish in the deconfined phase, but is sharply reduced. 



agrees with lattice calculations |T8[: The quenched chiral condensate does not disappear in 
the deconfined regime, but it is sharply reduced as one crosses the phase transition, cf. in 
particular Fig. ^. In view of the discussion of section |^, the fact that the smoothest choice of 
gauge field construction yields the most realistic physics is not unexpected. This appears to 
be the model of choice. As discussed above, the fact that it yields a value for the quenched 
chiral condensate in the confining phase which is slightly enhanced compared with lattice 
calculations has no direct phenomenological relevance. Alternatively, one could envisage 
ultimately using a minimal monopole model which interpolates between the smooth and 
random S cases such as to reproduce the lattice results even more closely. 



32 





-(gg)/GeV^ 


T/T, 


minim, monop. 
smooth S 


minim, monop. 
random S 


random monop. 
smooth S 


random monop. 
random S 





0.129 


0.108 


0.175 


0.100 


0.83 


0.094 


0.100 


0.157 


0.097 


1.1 


0.059 








1.4 


0.0033 








1.67 


0.0015 


0.078 


0.103 


0.203 



Table 1: Chiral condensate as a function of temperature for different gauge field constructions 
on a lattice of spatial extension 1.56 fm, obtained by extrapolating the bulk of the Dirac 
spectral density to zero eigenvalue. Since the inverses of the temperatures T = 1.1 and 
T = 1.4 Tc are not integer multiples of the lattice spacing a in the random surface model with 
curvature coefficient c = 0.24, the associated values of the chiral condensate were obtained 
by the following interpolation procedure |l[]. The inverse of T = 1.1 Tc corresponds to one 
lattice spacing for c = 0.14 and to two lattice spacings for c = 0.32; interpolating {qq) 
as a function of c defines the chiral condensate also at c = 0.24 for T = 1.1 Tc. The same 
procedure was also applied to T = 1.4 Tc, the inverse of which corresponds to one lattice 
spacing for c = 0.2 and to two lattice spacings for c = 0.407. 

6 Summary 

Three fundamental nonperturbative phenomena characterize the low-energy regime of the 
strong interaction: Confinement, the Ua{^) flavor anomaly, and the spontaneous breaking of 
chiral symmetry. While the {SU{2) color) random vortex world-surface model was shown to 
quantitatively describe the first two of the aforementioned phenomena in 0,0, the present 
work indicates that also the third phenomenon, i.e. the spontaneous breaking of chiral sym- 
metry, finds a viable description within this model. Specifically, a gauge field construction 
for vortex world-surface configurations was found which, coupled to quarks, leads to a behav- 
ior of the quenched chiral condensate as a function of temperature in qualitative agreement 
with the behavior found in lattice Yang-Mills theory. Quantitatively, the quenched chiral 
condensate (qq) differs from the lattice Yang-Mills result by a factor of the order of unity 
which presumably can be absorbed into the ultraviolet regularization scheme dependence 
of the chiral condensate, or, equivalently, into the definition of the current quark mass m. 
Only the product m(gg) is a renormalization group invariant, physical quantity which can be 
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tested e.g. using current algebra techniques. Stringent quantitative tests of low-energy quark 
physics via additional hadronic observables are deferred to future work. 

The main thrust of the present investigation lay in developing the technical tools necessary 
for such work. In order to incorporate quark physics into the vortex model, the Dirac operator 
in a vortex background had to be constructed. This involves, as a prerequisite, casting 
vortex configurations explicitly in terms of gauge fields. Having appropriately specified such 
vortex gauge fields via a Wu-Yang construction, a matrix representation of the Dirac operator 
in a truncated (finite element) quark basis was obtained which determines infrared quark 
propagation in the combined vortex-quark system. The evaluation of the quenched chiral 
condensate highlighted above represents a first exploratory application of these techniques. 

With the Dirac operator at hand, one can envisage also carrying out dynamical quark 
calculations within the vortex model, by reweighting the vortex ensemble with the Dirac 
operator determinant. This will substantially penalize Dirac eigenvalues of very small mag- 
nitude and thus presumably reduce the chiral condensate to the phenomenologically expected 
values around (230 MeV)'^. Moreover, to make contact with phenomenology, the random vor- 
tex surface model must still be extended to SU{3) color; also in the three-color case, lattice 



Yang- Mills measurements [36| indicate that vortex degrees of freedom dominate the infrared. 



confining sector of the strong interaction. 
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